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Abstract 

We determine the convergence radius m conv for the expansion in the cur- 
rent quark mass using the Dyson-Schwinger (DS) equation of QCD in the 
rainbow approximation. Within a Gaussian form for the gluon propagator 
DfivvP) ~ SpuX e A we fi n d that m conv increases with decreasing width 
A and increasing strength \ 2 ■ F° r those values of \ 2 an d A, which pro- 
vide the best known description of low energy hadronic phenomena, m conv 
lies around 2Aqcd, which is big enough, that the chiral expansion in the 
strange sector converges. Our analysis also explains the rather low value of 
mconv ~ 50 ... 80 MeV in the Nambu-Jona-Lasinio model, which as itself can 
be regarded as a special case of the rainbow DS models, where the gluon 
propagator is a constant in momentum space. 
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The chiral perturbation theory (xPT) has turned out to be a rather powerful technique 
for analyzing hadronic phenomena involving Goldstone bosons (tt,K) as well as nucleons at 
low energies []I|-£§. One crucial assumption in this approach is the smallness of the chiral 
symmetry breaking current mass mo compared with Aqcd ~ 200MeV. This is surely fulfilled 
in the u-d quark sector (m u = 5MeV, = 9MeV) but questionable in case of s-quarks 
(m s pa 150MeV) @]. It is therefore important to look at the convergence of the chiral 
perturbation expansion for values of quark masses in this region. There are mainly two 
features, which determine the convergence of this series: 

1. Non-analytic terms, such as mo 3 ' 2 or In mo, which arise from Goldstone boson (7r,K) 



2. The convergence of the power series c n mo n itself already on the mean field level, i.e. 
without any boson loops 0. 

In this paper we will focus completely on the second point. This problem has been studied 
originally in the framework of the chiral a model by Carruthers and Haymaker (CH) || and 
recently within the Nambu-Jona-Lasinio (NJL) model |7j by Hatsuda HIJ. The basic idea 
as well as the main results are practically the same in both cases Q Let us take the NJL 
model as an example and briefly review the main points. The convergence radius m conv is 
determined by the existence of a solution of the classical equation of motion (gap equation) 
for the constituent quark mass M, which reads in Euclidean space: 



where G denotes the strength of the four fermion point coupling and Kuv the ultraviolet 
cutoff. Eq. (H) is solved for M with as input, which gives a relation M = M(m ) (c.f. 
fig.!). Hereby G and Auv are fixed and adjusted to reproduce the experimental value of 



loops 0. 




(1) 



x One should note that the chiral a model and the NJL model are closely connected through the 
gradient expansion |1C]. 
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the pion decay constant f n as well as a constituent quark mass of M = 300MeV in the 
chiral limit (mo = 0), which can be assumed to be a physically reasonable value. Point (A) 
in figjl] corresponds to the vacuum with spontaneously broken chiral symmetry (mo = 0, 
M = 300MeV). This is the point around which the chiral expansion in mo is performed. As 
one can see the relevant region for the convergence of the series 

M(m„) = E c «™o n (2) 

n 

are the negative values of mo, i.e. the branch of the curve between the points (A) and 
(B). If m is smaller than —\m conv \ (point (£?), where |^ = 0) no solution of eq.(|l]) exists, 
that can be reached continuously from A. This means the breakdown of the series (@) for 
all values of m with |mo| > |m com) |. From fig.[3] we can read off m conv m 70MeV. It has 
been checked that this value does not depend crucially on the form of the regularization 
scheme (e.g. 0(3) and 0(4) sharp momentum cutoffs, proper time etc. ||11||). From this it 
follows that in case of the NJL model the strange quark mass m s « 150MeV lies far beyond 
m conv . Hatsuda has demonstrated the failure of the convergence of the series (0) explicitly 
for various observables such as the quark condensate (qq) or the nucleon sigma term 

EE- 

It is the aim of this paper to study this problem within the rainbow Dyson-Schwinger 



(DS) approach to QCD fl2| , |T3ll (for a comprehensive review c.f. [14|). This model has been 



extensively used for describing mesonic and vacuum properties [fq-Qq] and at least ten 



tatively also for the calculation of nucleon observables |T9H23|. The general form of the 
rainbow DS equation for the quark propagator S*(p) reads: 

H)£(P) = \ds 2 J ^[(iy»)(iS(k))(iD^(p - k)(i7„)] (3) 
where S(p) denotes the quark self energy defined by 

S- 1 (p)=^-m -S(p) (4) 

and D^lp — k) the gluon propagator. After continuation to Euclidean space time and 



using the Feynman gauge with the running coupling constant a(p 2 ) (where p denotes the 
Euclidean 4-momentum and p = ^/p ■ p) 

D^(p) = 5^D(p 2 ) = S^^^J- (5) 

g s 2 p 2 

as well as the decomposition of the self energy £ 

^(p) = ^A(p 2 )-l] + B(p 2 )-m (6) 
one obtains the coupled set of equations: 

r>W) - 1] = |*V t&^-d + *tf) (7) 

Btf) = + f / - ,)■) ^gg^ (8) 

As soon as a form for D(p 2 ) is specified the eqs. (0) and (|§|) define a certain model. It is 
easy to convince oneself that the gap equation of the NJL ([!]) appears as a special case of 
eqs. (0) and (0) if the gluon propagator is a constant in momentum space 

D NJL {P 2 ) = -7T-G (9) 



or equivalently 



«7VJL = ( 10 ) 



which gives a contact interaction in coordinate space with the momentum independent 
solutions A(p 2 ) = 1 and B(p 2 ) = M. From a principle point of view none of these models, 
i.e. none of the gluon propagators is preferred over the other, because none of them can be 
derived from QCD. On the other hand the general features of QCD are more likely better 
described with a running coupling a(p 2 ), which on the one side grows in the infrared region 
and therefore gives rise to "confined" quarks and on the other side shows asymptotic freedom 
in the ultraviolet region 

2 p^oo __dm 

In (*7 W 
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(with d 



12 



33-2N f 'I 



rather than with a form like olnjl{p 2 ) ~ P 2 , which gets small in the IR 
region but grows for high p 2 and has to be cut off at Aqcd- In refs. |]T5| , p~7|| a superposition of 



a Gaussian form with strength x 2 and width A (dominating at small p 2 ) and the standard 
asymptotic form (dominating at high p 2 ) 



a(p 2 



3nx 2 ( p 



4 



A 2 + 



dn 



(12) 



ln(r + p 2 /A 2 ) 

with A = Aqcd = 190MeV, x = 1.14GeV, A = 0.002GeV 2 and r = 3.0 has been shown to 
provide a reasonable description of low energy hadronic phenomena. Furthermore it turned 
out that the strength of a(p 2 ) in the IR domain is large enough that the model shows 
spontaneously broken chiral symmetry, i.e. the DS equation (0H) nas a nontrivial solution 
B{p 2 ) ^ if m = 0, as well as "confinement", in the sense that the quark propagator S(p) 
does not have a pole at timelike p 2 . 

For our study of the convergence of the chiral expansion in this model we will first 
consider a pure Gaussian ansatz 



a{p 2 



3tt X 2 ( p 2 



-P 2 /A 



(13) 



4 \A 2 / 

and vary the strength x 2 as well as the width A. For the determination of m conv we have 
again to look at negative m . In order to do so we write the set of the DS equations (^,||) 
in the form: 



4 r c 

P 2 [A(P) - 1] = r / 
3 Jo 



; A{q 2 ) (pq) 
q 2 A 2 (q) + B 2 (q) 



£a{p, q) 



B(p) = B(0) + - / dqq 
3 Jo 



B(q) 



q 2 A 2 (q) + B 2 (q) 
where the integral kernels Ka and Kb axe given by: 



[K. B {p,q)-K B $M 



(14) 
(15) 



K>a{p,q) = J dVt 



q(p - qf 
(p - q) 2 



^ Pq ) = 4A 6 



^ \a) +i 2 {-a 



(16) 



and 



K B {p,q)= dn 



g(p-qY 
(p - q) 2 



(17) 



respectively. Hereby I n denotes the Bessel function of order n. The set of integral equations 
(|Tj]) and ([0|) is formally independent of mo but instead depends on the initial value B(0). It 
can be uniquely solved for any given -B(O) and the corresponding mo can then be extracted 
from the solution B(p) at high p 

lim B(p) = mn (18) 

p— *oo 

This renders a relationship between -B(O) and m , i.e. m = m [-B(0)], analogous to the one 
displayed in fig.|l[ The convergence radius m conv is then determined by the minimum of this 
curve, i.e. the point: 

dm 



dB(0) 



(19) 



The set of eqs.(|14]) and (|TJ) is solved by a selfconsistent procedure in the interval 
[0, 1000A] with logarithmic grid points. 

1. Let us first look at solutions with various widths within the Gaussian parameterization 
( P^D keeping x 2 fixed at x — 8. OA. The forms of the corresponding running coupling 
constants a(p) are displayed in fig.^|, where for comparison we have included the one of 
the NJL eq.fllOD in addition. In fig.||] we compare the solutions B(p) for a fixed value of 
-B(O) = 7.75, whereas fig. |] shows the dependence m = m [£>(0)]. It is interesting to 
look at the limes A — ► 0, which means, in fact, that the gluon propagator approaches 
a 5-function in momentum space (dotted line in fig.f|): 

D{j> - q) = ( 27 r) 4 Ai^(4) (p _ q) (20) 
16 g s 2 



In this case we obtain an algebraic relation between m and B(0) from eq.([8|): 

y2 

m ° = m ~ W) (21) 

From figs.^l and ^ we can clearly see, that m conv rises with decreasing A. For small 
values of A (e.g. A = 0.2 A 2 , full line in fig||) one faces the problem that the self 
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consistent iteration procedure for the solution of the system fll4|) and fll5|) does not 
converge but ends up switching between two or more configurations, if the value of 
-B(O) deviates too much from the one in the chiral limit (mo = 0). This is a typical 
feature encountered in many nonlinear systems and does not mean that there exists 
no solution at all, but only that the simple selfconsistent procedure is not able to find 
it. It might be possible that more elaborate methods are successful in this case. On 
the other side we are clearly able to give at least a minimum value for m conv from 
the corresponding curve in fig.f|, which is sufficient for our analysis. Furthermore we 
recognize that for those small A the curve, as far as it can be calculated from the 
self consistent procedure, lies very close to the one obtained with a 5-function gluon 
propagator (A = 0, dotted line in fig.|). 

2. For A = it is clear from eq.(^l|) that 



{m [B(0)]} xi > {m [B(0)}} 



\2 



(22) 



if xi < X2- We have checked that this relation holds generally if the value of A. From 
this we conclude that m conv increases with increasing strength x 2 - 

3. Finally we have considered also gluon propagators of the form (0), which have in 
addition to the Gaussian form in the IR region have the logarithmic tail at high p 2 
leading to asymptotic freedom. In this case the momentum integrals in eqs. ([?]),(§), 
( |14D and ([15]) run from to the "renormalization point" fi 2 and the m is considered as 
the running mass m (/i). One can convince oneself, that the DS equations eqs.(|7]) and 
(^D or (|14D and (|T5| ) are consistent with the renormalization group results p4|- p6|JT3[] : 



B(p 2 ) 



2\ 4 



m (n) In 



A 2 



2\-d 



111 



p_ 

A 2 



(23) 



and 
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2 



m [fi) = m (/i) -^j (24) 

The general behavior of the relationship mo(/i) = mo(/x)[-B(0)] is the same as without 
this asymptotic tail. With the parameter set from refs. [JH| , |IT| mentioned above we 
obtain a convergence radius of at least m conv (fi = 1000A) = 1.6A. Due to eq.(|24]) this 
corresponds at a typical hadronic scale of /i ~ lGeV to a value of m conv (fi ~ lGeV) ~ 
3A, which lies clearly above the strange quark mass. 

Our results can be summarized as follow: 

1. In the framework of the rainbow DS approach, where the form of the gluon propagator 
D(p 2 ) is taken as input, the convergence radius of the chiral perturbation expansion 
mconv turns out to be crucially dependent on the infrared behavior of the gluon prop- 
agator. 

2. m conv gets larger if in the IR domain the overall strength of the gluon propagator 
increases or its width decreases. 

3. Using a parameterization D(p 2 ) which provides a good description of hadronic phe- 
nomena at low energies we obtain a value for m conv which is clearly larger than the 
strange quark mass and therefore the chiral expansion in the strange quark sector 
converges. 

4. Especially we are now able to explain the rather low value of m conv and therefore the 
poor convergence of the chiral expansion in the strange quark sector in case of the NJL 
or the chiral a model, where the corresponding gluon propagator is a "small" constant 
in momentum space. Our results indicate that this small value of m conv ~ 50 . . . 80MeV 
in these models arises more likely due to the special form of a contact interaction than 
it is a general feature of QCD and the early breakdown of the chiral expansion does 
not occur in a theory with infrared slavery. 
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FIGURES 

FIG. 1. Determination of m conv in the NJL with proper-time regularization, Ajjy = 634MeV, 
G = 30A UV 2 . 



FIG. 2. Comparison between running coupling constants a(p): Gaussian form (eq.(|13|) with 
X = 8A and A = 200MeV) for various values of A and the NJL (eq.(|lO|) with A uv = 634MeV, 
G = 30A UV 2 ). 



FIG. 3. The solution B(p) of the DS equations ( |14l ) and (15) for a Gaussian running coupling 
a(p) with various widths A ( X = 8A, B(0) = 7.75A, A = 200MeV 

)■ 



FIG. 4. The dependence tjiq = mo[B{0)] obtained from the solutions of the DS equation (13) 



using a Gaussian running coupling compared with the <5-fct. limes A = (eq.(^)) and the NJL. 
The parameters are the same as in fig.0. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405289vl 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405289vl 
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This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405289vl 
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This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405289vl 



